MOMENTUM OPERATORS IN TWO INTERVALS: 
SPECTRA AND PHASE TRANSITION 



PALLE E. T. JORGENSEN, STEEN PEDERSEN, AND FENG TIAN 

Abstract. We study the momentum operator defined on the disjoint union 
of two intervals. Even in one dimension, the question of two non-empty open 
and non-overlapping intervals has not been worked out in a way that extends 
> the cases of a single interval and gives a list of the selfadjoint extensions. 

Starting with zero boundary conditions at the four endpoints, we characterize 
C'~') the selfadjoint extensions and undertake a systematic and complete study of 

the spectral theory of the selfadjoint extensions. In an application of our 
extension theory to harmonic analysis, we offer a new family of spectral pairs. 
Compared to earlier studies, it yields a more direct link between spectrum and 
geometry. 
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1. Introduction 

An early paper of John von Neumann [vN32] deals with unbounded linear op- 
erators in Hilbert space, and their spectral theory. By now the theory is part of 
functional analysis, and the problems entail the study of "deficiency indices;" now 
important tools in mathematical physics; see e.g., [Kre49]. 

Indeed, von Neumann's initial motivation came directly from quantum mechan- 
ics: measurements of observables, prepared in states. These indices have now be- 
come the accepted framework in spectral theory, as well as in quantum mechanics. 

In quantum theory, states are realized by vectors (of norm one) in Hilbert space 
J^ 5 , and the quantum observables correspond to selfadjoint operators in ,#C (the 
state space.) 

Our system of operators in two intervals has a number of applications not con- 
tained in other studies. For example, when modeling quantum states (Section 4), 
we show how a two-component quantum system offer a binary model for the dynam- 
ics of wave functions. Indeed the parameters for our family of selfadjoint operators 
allow a direct computation of transition probabilities. 

Here we study an instance of the index/spectral theory question, and in compu- 
tations we use the spectral theorem as it is applied to selfadjoint operators. However 
the deficiency indices are precisely the obstruction to selfadjointness for unbounded 
Hermitian operators. 

This is illustrated with examples, naturally occurring, and quite simple, from 
quantum theory: There one studies wave-functions with confinement to a finite 
region (particles in a box), e.g., confinement in a finite interval. In these cases, the 
physical operators are not selfadjoint, but only Hermitian. 

But the unbounded linear operators will still have dense domain in the physical 
Hilbert space (of states). Therefore, an assignment of boundary conditions will 
amount to a subtle extension of the initial dense operator domain. 

In several variables, one is then looking for commuting selfadjoint extensions; see 
[Fug74, J0r82, PW01, Ped04, JPOO, JP99, JP98, DJ07] for a study of this problem, 
in several applications. 

For these problems, and for applications in the theory of linear partial differen- 
tial operators, one is led to a more axiomatic consideration of boundary conditions. 
Indeed, J. von Neumann, in his paper [vN32], identified such conditions in the 
Hilbert space theoretic framework. For a single Hermitian operator, they take the 
form of deficiency spaces, and deficiency indices, see details below. In one dimen- 
sion, and for a finite interval, say I, a popular model example has its deficiency 
indices come out (1,1). Hence, by von Neumann's theorem, this yields a correspond- 
ing one-parameter family of selfadjoint extensions. In the case of one interval, they 
can readily be computed, and it is possible to read off the spectrum of each of the 
selfadjoint extensions. Each one is the generator of a unitary one-parameter group 
U(t) acting in Jif. And each of these unitary groups of operators in turn translates 
the wave-functions in the interior of the interval, and then assigns a phase-factor at 
the endpoints of /. The corresponding dynamics of the wave-functions turns into 
a fixed phase-shift, i.e., a matching with a phase-factor of incoming and outgoing 
waves. 

Moreover, for each selfadjoint extension, the corresponding spectrum is a Fourier 
spectrum for the Hilbert space — L 2 (I) where I is the given interval. By "Fourier 
spectrum" we mean that the spectrum (of the selfadjoint extension) is discrete of 
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uniform multiplicity one, and that the eigenfunctions are complex exponentials. We 
say that the pair, set and spectrum, is a spectral pair. See [JP98] for the theory of 
spectral pairs. 

Indeed, for more general models with wave equations, and in higher dimensions, 
there is a scattering theory model which accounts for scattering of waves at obsta- 
cles, see for example [LP 68]. Yet, still staying with the translation model, but going 
from one dimension to M. d , Fuglede [Fug74] asked for an analogue of the translation 
model in one dimension, but now for bounded open domains in R d . 

Nonetheless, even for one dimension, the question of two non-empty open and 
non-overlapping intervals has not been worked out in a way that extends the cases 
of a single interval and gives a list of the selfadjoint extensions. 



1.1. Two Intervals. This paper is devoted to the two-interval case. Now, starting 
with zero boundary conditions at the four endpoints, we show that the deficiency 
indices are then (2,2), see [Kre49]. Moreover, we offer a characterization and pa- 
rameterization of all the selfadjoint extensions. For each one, we compute the 
spectrum, and we write down the corresponding unitary one-parameter groups. 

We show that, by contrast to the one-interval case, that some of the selfadjoint 
extension may have points in the spectrum of multiplicity two. 

Our analysis lies at the crossroads of operator theory, spectral geometry, and 
harmonic analysis. For the two-interval case, we find the selfadjoint extension 
operators Pb, and their spectra. Each operator Pb, indexed by the group U(2) of 
unitary 2 by 2 matrices, has pure eigenvalue spectrum {X n (B)}; the multiplicity 
can be 1 or 2. In fact, we have the eigenvalues represented as functions of both B 
and of the choice of endpoints in I 2 ■ 

We study all possible configurations of the two intervals. It turn out that only the 
ratio r of their lengths is important. The nature of this number has subtle spectral 
theoretic implications; e.g., the distinction between the cases when r is rational vs 
irrational is important. When it is irrational, the points in the spectrum of Pb 
form a generator for pseudo-random numbers [BAAS11]; see Section 3 for details. 
For the general case, we display the points {A„} in the spectrum of Pb as functions 
of 5 variables, 4 for B in 17(2), and one for r. 

In Section 5, we determine completely the cases when we the two intervals, 
together with spectrum of Pb, yield spectral pairs (SP). It turns out to be a small 
subset of the total parameter space. In the SP-case, we show that the set A of points 
{A„} fall in two simple cases: Modulo the angular part of B, A is either is the union 
of one or two lattice co-sets. For simplicity suppose the intervals are [0, 1] U [a,/3] 
where 1 < a < (3. In Theorems 5.3 and 5.8, we give a complete analysis of the SP 
cases. Modulo angular variables, this is merely a "small" discrete subfamily: Firstly, 
SP forces either (3/(1 + ft — a) to be an integer or j3 — a = 1; and secondly, we get 
the eigenvalue lists as follows: Let w be the radial variable in B. If (3/(1 + (3 — a) 
is an integer then the points \ n form a lattice, and if /? — a = 1 the points A n (as a 
function of w) contains the branches of the arc-cos function. In summary, the SP 
cases form a singular family, representing a broken symmetry, or phase transition; 
a discrete family within a much bigger continuum. For comparison, note that in 
[Rob71], dealing with second order operators in a single interval, the author finds 
an analogous spectral analysis of phase transitions. 
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1.2. Unbounded Operators. We recall the following fundamental result of von 
Neumann on extensions of Hcrmitian operators. 

Lemma 1.1 (see e.g. [DS88]). Let L be a closed Hermitian operator with dense 
domain f^o in a Hilbert space. Set 

@± = {4>± e dom(L*) | L*i>± = ±#±} 
tf(L) = {U : ®+ -> 9- \U*U = P@ + ,UU* = P®_] (1.1) 
where P@ ± denote the respective projections. Set 

g(L) ={S\LCS,S* = S}. 

Then there is a bijective correspondence between ^(L) and <£(L), given as follows: 
If U G C &{L), and let Ljj be the restriction of L* to 

{ipo + f+ + Uf+ | cpo G %, f+ G (1.2) 

Then Ljj G £"{L), and conversely every S G S(L) has the form Ly for some 
U G With S G tafce 

U-iS-UXS + il)- 1 ^ (1.3) 

and note t/iat 

(1) C7 G ^(i), and 

(2) S = L V . 

Vectors f in dom(L*) admit a unique decomposition f = ipo + /_)_ + /_ where 
ipo G dom(L), and f± G For £/ie boundary-form B(-, ■), we fta?;e 

»B(/,/) = (L*f,f)-(f,L*f) 

= ||/ + f-||/-|| 2 . 

1.3. Summary. We undertake a systematic study of interconnections between ge- 
ometry and spectrum for a family of selfadjoint operator extensions indexed by 
two thing: by (i) the prescribed configuration of the two intervals, and by (ii) the 
von Neumann parameters. This turns out to be subtle, and we show in detail 
how variations in both (i) and (ii) translate into explicit spectral properties for the 
extension operators. Indeed, for each choice in (i), i.e., relative length of the two 
intervals, we have a Hermitian operator with deficiency indices (2,2). Our main 
theme is spectral theory of the corresponding family of (2, 2)-selfadjoint extension 
operators. Accidentally, this more general problem has implications for a question 
from [LabOl] on spectral pairs. But our study of the spectral theory dictated by 
(i) and (ii) goes far beyond its implications for spectral pairs. 

Another advantage of our approach via a definite 6-parameter family of selfad- 
joint operators is that the verification for completeness of the associated function- 
systems becomes easy. We first prove that all these operators have pure point- 
spectrum, see Theorem 2.4. In each verification, one may then take advantage of 
the spectral theorem applied to the particular selfadjoint operators under study. 
So this applies whether or not one has a spectral pair. Otherwise the completeness 
property is typically subtle in spectral geometry and in physics applications; see 
e.g., [DS88, Gil72, KM06, LabOl, tab02, LP68, Rob71]. 

One of our applications will be to the characterization of spectral pairs in the 
two-interval case(s). While this question was addressed in the paper [LabOl], it 
was attacked there with the use of geometric tools; specifically by listing the ways 
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two intervals can tile the real line under translations. We obtain the same results 
regarding spectral pairs, but now with the use of operator /spectral theory. Even 
though most of our computations are for interval examples, we believe they are 
of general interest, and that they offer intriguing extensions to more structured 
models. 

1.4. Prior Literature. There are related investigations in the literature on spec- 
trum and deficiency indices. For the case of indices (1, 1), see for example [ST10]. 
For a study of odd-order operators, see [BH08]. Operators of even order in a single 
interval are studied in [Oro05]. The paper [BV05] studies matching interface con- 
ditions in connection with deficiency indices (m, m). Dirac operators are studied in 
[Sak97]. For the theory of selfadjoint extensions operators, and their spectra, see 
[Smu74, Gil72], for the theory; and [Naz08, VGT08, Vas07, Sad06, Mik04, Min04] 
for recent papers with applications. For applications to other problems in physics, 
see e.g., [PR76, Bar49]. 

2. Momentum Operators 

By momentum operator P we mean the generator for the group of translations in 
L 2 (— oo, oo). There are several reasons for taking a closer look at restrictions of the 
operator P. In our analysis, we study spectral theory determined by a choice of two 
intervals (details below.) Our motivation derives from quantum theory, and from 
the study of spectral pairs in geometric analysis; see e.g., [DJ07], [Fug74], [JP99], 
[LabOl], and [PW01]. In our model, we examine how the spectral theory depends 
on both variations in the choice of the two intervals, as well as on variations in the 
von Neumann parameters. 

Granted that in many applications, one is faced with vastly more complicated 
data and operators; nonetheless, it is often the case that the more subtle situations 
will be unitarily equivalent to a suitable model involving P. This is reflected for 
example in the conclusion of the Stone- von Neumann uniqueness theorem: The 
Weyl relations for quantum systems with a finite number of degree of freedom are 
unitarily equivalent to the standard model with momentum and position operators 
P and Q. 

2.1. Applications. In Lax-Phillips scattering [LP68] theory for the acoustic wave 
equation, one obtains the solution to a particular wave equation represented by 
a unitary one-parameter group acting on Hilbert space. Wave equations preserve 
energy; hence the Hilbert space. 

With scattering against a finite obstacle in some number d of dimensions, one 
then gets two closed subspaces corresponding to incoming and outgoing waves, 
respectively, here referring to the whole space, in this case K d . The Lax-Phillips 
scattering operator is derived from the two. Again, in the two subspaces, the 
respective unitary one-parameter groups are unitary equivalent to a P-model. By 
taking into account uniform multiplicity, one arrives at the operator P. 

For other uses of P and its restrictions, see [J0r81] and [Rob71]. The application 
in [Rob71] is to quantum statistical mechanics. 

In geometric analysis, one studies second order Hermitian PDOs on a finite 
region in K d , and selfadjoint boundary condition, for example Dirichlet conditions, 
or Neumann conditions, see also [LP68], [DS88] and [Rob71]. 
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In quantum systems, the initial interest is in the Schrodinger equation for the 
dynamics of quantum states. And once again an essential part of the problem is 
turned into the study of a unitary one-parameter group in Hilbert space govern- 
ing the solution to the Schrodinger equation: We are faced with determining the 
spectrum of the selfadjoint generator. Even if there may not be a direct unitary 
equivalence, there will often be a part of the dynamics, e.g., bound-states, local- 
ization, and quantum tunneling, where the two-interval model offers insight. In 
[Rob71] for example, the operator is second order, and it is used in the study of 
confinement quantum states (wave functions) in a single interval. There is a con- 
nection between this, and the details below. Below we have a first-order operator, 
and two intervals. But in both instances, the models are represented through the 
study of deficiency indices (2, 2). In both cases, one is able to determine an eigen- 
value spectrum (i.e., bound states) from a system of curves with intersection, but 
the systems vary from one application to the other. For details, in [Rob71], see 
Figures 1 through 8. The analogous system of curve intersections of our present 
two-interval models are given in Examples 3.21 and 5.15 below. 

In studies of quantum scales, one is faced with quantum tunneling, see e. g., 
[Rob71], a process that cannot be directly perceived at a macroscopic scale. Quan- 
tum particles (or rather waves) travel between potential barriers. With a small 
probability, wave functions will crossing barriers. See Corollary 4.1 for a discus- 
sion. The reason is duality for quantum states; i.e., states having simultaneously 
properties of waves and particles. This leads to computation of probability density 
for particle's position, thus describing the probability that the particle is at any 
given place. In the limit of large barriers, the probability of tunneling decreases for 
taller and wider barriers. The simplest and best understood tunneling-barriers are 
rectangular barriers, offering approximate solutions. 

2.2. The boundary form, spectrum, and the group U(2). Since the problem 
is essentially invariant under affine transformations we may assume the two intervals 
are 1\ := [0, 1] and I 2 := [cx,0\. L 2 {I\ U I 2 ) is a Hilbert space with respect to the 
inner product 

(/Is) ■= f fg+ [ fg- (2.1) 
Jii J12 

The maximal momentum operator is 

^=44 (2-2) 

i2tt dt y ' 

with domain 3{P) equal to the set of absolutely continuous functions on ii U /2. 
The boundary form associated with P is the form 

B(f,g):=(Pf\g)-{f\Pg) (2.3) 

on !&(P). Clearly, 



B(/, 9) = /(1)<7(1) - /(0).g(0) + f{P)g{[3) - f(a)g(a). (2.4) 

For / € 9(P), let Pl (f) := (/(l),/(/3)) and p 2 (f) := (/(()), /(a)) . Then 

B(/,. 9 ) = (pi(/) I pi{ 9 )) - (p 2 (f) I P2{g)) ■ (2.5) 

Hence (C 2 , pi, p 2 ) is a boundary triple for P. The set of selfadjoint restrictions of 
P is parametrized by the set of unitary 2x2 matrices, see e.g., [dO09]. Explicitly, 



MOMENTUM OPERATORS IN TWO INTERVALS 



7 



any unitary 2x2 matrix B determines a selfadjoint restriction Pb of P by setting 

9 (Pb) :={fe@ (P) | B Pl (f) = p 2 (f)} . (2.6) 

Conversely, every selfadjoint restriction of P is obtained in this manner. For com- 
pleteness, we begin by working out the details connecting the boundary form for- 
mulation to the von Neumann deficiency space approach in our setting. 

Definition 2.1. For I\ and I2, consider the following two copies of the two- 
dimensional Hilbert space C 2 : 

P2 {9{P)) = m L = U ) |/g#(p) 



Pl (9{p)) = a R = U {ft! ) / = '/'{P ) 



and 

n . (Oil p^^ — — ) I 

fW) 

where p\ and ^2 are the respective boundary-restrictions. 

For / G @(P), let / = (^0 + /+ + /- be the decomposition in (1.2), i.e., ^ € ^0, 
/± G $t±. Note that, for / <E @(P), we have 

B(/,/) = H/flll^ - ||/ L || 2 ^ = II/+H 2 - ||/_|| 2 . (2.7) 

Set 

% = {B : — \ S) R , isometric}, 
i.e., a copy of the matrix-group of all unitary 2x2 matrices. 

Proposition 2.2. In the two-interval example, there is a natural isomorphism 
c &b — C &{L) where ^(L) is defined as in Lemma 1.1. 

Proof. Note that by Lemma 1.1, each of the selfadjoint extensions from the minimal 
domain ^0 (zero-boundary conditions) is determined by some U € ^(L) via /_ = 
Uf+, /+ G @+. Specializing (2.7) to the domain %j in (1.2), i.e., / = tp +f + +Uf+, 
(fio G %, f+ G 9 + , we get 

B(/,/) - H/flJI 2 - H/lII 2 = II/+H 2 - \\Uf+\\ 2 = 
for all f £ 3>u- It follows that 

h ^ /*, (2.8) 
defined for / G &u , defines a unique element Btfj, and that 

tf(L) 3U^ B e% (2.9) 

is an isomorphism. ■ 

Since the two deficiency spaces @± are easy to compute in each of the two- 
interval examples, we can write down two isomorphisms 2$+ -4- and ^?r; 
note that all four spaces are two-dimensional. With this, we get the following 
representation of the isomorphism in Proposition 2.2: 

Corollary 2.3 (Matrix realization). Let U, b±, and B be the operators in Propo- 
sition 2.2, then B : 3$i, — > £3 R is B = b-Ub^ 1 ; see Figure 2.1. 
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L 

B=b-Ub~ 
R 



Figure 2.1. 

With (2.9) in Proposition 2.2, we now parametrize the collection of unitary 2x2 
matrices by 

g _ ( w e((p) —yl — w 2 e(9 — tp) 
\ Vl — w 2 e{tp) w e(8 — <fi) 

where < w < 1, <f>, ip, 9 € R, and 

e(ar) := e l27ra (2.11) 
Some of our results about the spectrum are summarized in 

Theorem 2.4. If Pb = Pw.4>,ip,e is the selfadjoint momentum operator on [0, 1] U 
[ot, 0\ associated with the unitary matrix (2.10) via (2.6), then Pb has pure point 
spectrum. In fact, the spectrum of Pb is the set of all real solutions A to the equation 

(e {(j> + A) - w) e (9 - <f> + (/? - a) A) = w e (cf> + A) - 1. (2.12) 

In particular, the spectrum is independent of the parameter ip and does not 
depend on the location of the interval [a, f3) only on its length. 

Remark 2.5. Our analysis also works when a — 1. In this case one has to be careful 
with the interpretation of the function values at the point of intersection of [0, 1] 
and [a, 0\. One way is the interpretations 

/(l) = = lim/W, and f(a) = /(1+) = lim/(t). 

We leave the details to the interested reader. 

Remark 2.6. For applications, the role of the four parameters in (2.10) is taken up 
in Section 4 below. Here we merely add a note about the geometry entailed by the 
two extreme cases (i) w = 0, and (ii) w = 1. In the first case, the matrix B is 
off-diagonal, with zeros in the two diagonal slots, while for (ii) it is diagonal. As 
a result, for (ii), we are merely dealing with the orthogonal direct sum of separate 
boundary value problems for the two intervals I\ and 1% treated in isolation; so two 
separate index (1, 1) problems. By contrast, for (i), the two intervals get exchanged 
each time a local translation reaches a boundary point. As noted in Section 4 below, 
the dynamics in the intermediate case < w < 1, may be computed with the use 
of a binary random-walk model. 

3. Spectral Theory 

This section is about a family of selfadjoint operators Pb arising as restrictions 
to the union of two intervals of the maximal momentum operator. We explore how 
these selfadjoint operators Pb depend both on the prescribed boundary matrix B 
from eq. (2.10), as well as the two intervals, i.e., on the choice of the two numbers 
a and /?, see eq. (2.1). So the operators and their spectra depend on six parameters 
in all, and this dependency is charted in complete detail below. 
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Note that the difference j3-a is the length of the second interval. Further recall 
that the unitary 2 by 2 matrix B has four parameters, see (2.10); and that its radial 
parameter w lies in the closed interval from to 1. We say that the endpoints 
constitute the extreme cases, and the open interval < w < 1, the generic case. 
The role of w is analyzed further in Section 4 which deals with scattering theory 
for the unitary one-parameter groups generated by each of the operators Pg. 

In the present section, we show that each operator Pg has pure eigenvalue spec- 
trum {A„(P)} with the eigenvalues depending on all parameters, i.e., on both a 
and f3, and on the matrix B. This dependence is detailed below: In Section 3.1 
we deal with two extreme cases for B; Section 3.2 the generic case; and in Section 
3.3 we deal with the dichotomy for the value of the length of the second interval, 
rational vs irrational. 

We saw in Corollary 2.3 that every self adjoint restriction of the maximal mo- 
mentum operator P has the form Pb for some 2x2 matrix B as in (2.10). We begin 
our investigation of how the spectrum of Pb depends on the parameters 1 < a < /?, 
4>, tp, 9, and < w < 1 by stating two simple lemmas. 

Fix 1 < a < p. Let h := [0,1], and I 2 := [a,0\. Let e x (y) := e(xy). Any 
complex number A is an eigenvalue of the maximal momentum operator P in (2.2). 
For each complex number A the formal solutions to Pf = Xf are 

e x ( x ) '■= {a>Xli(%) + bxi 2 {x))e x {x), a,beC. (3.1) 

Lemma 3.1. The maximal operator P has spectrum equal to the complex plane 
and each point in the spectrum is an eigenvalue with multiplicity two. 

Proof. Since e<f ' b) is in L 2 (h U I 2 ), it follows from Pe^' b) = Xe { ^ b) and (3.1) that 
each A € C is an eigenvalue of P with multiplicity two. ■ 

Lemma 3.2. The spectrum of any selfadjoint restriction P of P equals the set of 
eigenvalues of P and each eigenvalue has multiplicity one or two. 

Proof. By Lemma 3.1 the spectrum of P is the A for which e^' h ^ is in the domain of 
P for some choice of constants a, b. The multiplicity is determined by the dimension 
of this set of constants. 

Suppose Pb is the selfadjoint restriction of P determined by the matrix (2.10) 
via (2.6). When the parameters for B are fixed, the expression on the right-hand 
side in eq. (3.1) is a function of the x-variable. The constants a and b in front of the 
indicator functions depend on P, and, in turn, they determine the eigenfunctions, 
and therefore the spectrum, of the associated selfadjoint operator Pg . 

Suppose Pb is the selfadjoint restriction of P determined by the matrix (2.10) 
via (2.6). Then the boundary condition Bpi(f) = p 2 {f) states 

/ we(4 >) -VT^e(6-4,) \f f(l) \ ( /(0) \ 

Consequently, A is an eigenvalue for Pg iff there is are complex numbers a, b such 
that / = e^ 6 ' satisfies the boundary condition (3.2), i.e., 

/ we(4 >) -VI - w 2 e{6 - ip) \ f oe(A) \ _f a 
\ \/l-w 2 e(V) we{9-4>) ) \ b e(X(3) ) ~ \ b e(Aa) 

Since Pg is selfadjoint any eigenvalue A must be a real number. ■ 
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Below we show that the real solutions A to (3.3) is a discrete set having no other 
accumulation points than plus/minus infinity. 

3.1. Boundary Cases. Our analysis depends on the parameter w. We begin by 
considering the extreme cases w — 1 and w = 0. 

Theorem 3.3 (w = 1). Fix 1 < a < (3. Let I\ = [0,1], and I 2 ■= [a, 0\. Suppose 
Pb be the selfadjoint restriction associated with (2.10) via (2.6) and w — 1. Then 
Pb has pure point spectrum Ai U A2, where 

Ai = -0 + Z (3.4) 

A 2 = I + - Z 3.5 

p — a p — a 

and corresponding eigenfunctions are 

4>+n = X[o.i]e0+„, and e ^J = X[a,p\e*±n.- (3-6) 

Consequently, each eigenvalue has multiplicity one, except the eigenvalues, if any, 
in Ai H A2 have multiplicity two. 

Proof. In this case the boundary condition (3.3) takes the form 

e(<f>) \ / oe(A) \ _ / a \ , , 

e(e-(f>) )\be(\p) )~\be{Xa) J' { ^' ) 

That is 

ae(<fi + X) = a 
be (9- + /3A) =be[aX) . 
The stated conclusions are now immediate. ■ 

Remark 3.4. It is easy to see that (3.4) and (3.5) are the solutions to (2.12) when 
w = 1. 

We present a few examples in order to illustrate some of the possibilities. 

Example 3.5. If f3 — a = 1 and 9 = 2<j> the spectrum equals —0 + Z and has 
uniform multiplicity equal to two. 

Example 3.6. If <f> = 0, 9 is rational and j3 — a is irrational, then the spectrum has 
uniform multiplicity equal to one. 

Example 3.7. If <j) = 9 = and ft — a = 2, then the spectrum is ^Z, the eigenvalues 
in Z have multiplicity two and the non-integer eigenvalues have multiplicity one. 

Example 3.8. If = ^ = and /3 — a is irrational, then has multiplicity two 
and the other eigenvalues have multiplicity one. 

In the other extreme case we have 

Theorem 3.9 (w — 0). Fix 1 < a < (3. Let I 1 = [0,1], and I 2 := [a, j3]. Suppose 
Pb be the selfadjoint restriction associated with (2.10) via (2.6) and w = 0. Then 
Pb has pure point spectrum 
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with uniform multiplicity equal to one and an eigenfunction corresponding to the 
eigenvalue jq-g~ , neZ, is 

1 (i +n )p-(l-a)0\ \ , s 

2^ + U 1+V* j^ + ^j £ te- (3 ' 9) 

Proof. In this case the boundary condition (3.3) takes the form 

-e(e-i>) \ ( ae(X) \ _ / a \ , , 

e(^) )\ be(Xp) J ~ \ be{Xa) J ' ^' iUj 

Writing (3.10) as two equations we get 

-b e{9 - tP + X/3) = a (3.11) 

ae(V> + A) = be(Xa). (3.12) 

The first equation shows that \a\ = \b\ . In particular, we may set 6=1 and 
a = —e{9 — ip + A/3). Then the eigenvalues are determined by the second equation 
which now states 

-e(0 + X + pX) =e(Aa). 
Hence the eigenvalues A are determined by 

6 + A(l + - a) e - + Z. 

The stated conclusions are now immediate. ■ 

Remark 3.10. It is easy to see that (3.8) are the solutions to (2.12) when w = 0. 

3.2. Generic Cases. We now consider the non-extreme cases < w < 1. We begin 
by showing that the eigenvalues are the solutions to (2.12), thereby completing the 
proof of Theorem 2.4. 

Lemma 3.11. Fix 1 < a < p. Let I\ = [0, 1], and I2 := [cc, /?]. Suppose Pb is the 
selfadjoint restriction associated with (2.10) via (2.6) and < W < 1. TTien a point 
X inR is an eigenvalue of Pb if and only it is a solution to the equation 

e(e-t + (P-a)X) = W ;.^-l (3.13) 
e (0 + A) — w 

in terms of the argument this is 

6-<p+(P-a)X= — arctan ( (1 ; ^) ^ (2,(0 + A)) \ ^ 

^ VP ; 2tt V 2w - + w 2 )cos(2tt(0 + A))/ v ; 
When X is a real solution to (3.13) then any corresponding eigenfunction is a mul- 
tiple of 

4 ( x ) = ( a Xh(x) + Xh( x )) e x( x )i ( 3 - 15 ) 

where a is determined by (3.18) below. In particular, the spectrum has uniform 
multiplicity equal to one. 

Proof. The eigenvalue equation (3.3) is 



a w e(4> + A) - by/ 1 - w 2 e{6 -ijj + pX)=a (3.16) 
a\/l - w 2 e(tp + A) + b w e(0 - <f> + PX) = b e(aX). (3.17) 

Equation (3.16) shows a = ^=> 6 = 0. So, we may set 6 = 1, which gives 
(3.15). 
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Solving both (3.16) and (3.17) for a, shows a equals 

Vl - w 2 e{6-j) + p\) _ l-we(9-<t> + (j3- a)X) 
w e (</> + A) - 1 Vl - w 2 ety + A - aX) 

Cross-multiplying gives 

(1 - w 2 ) e (0 + A + {fi - a) A) 
= {l-we{e-4>+{P- a)X)) (w e (<f> + X) — 1) . 
Expanding the right hand side we see 
(l — w 2 ) e (9 + X + ((3 — a) X) 

= -1 - w 2 e {9 + X + Q3 - a)X) + w (e (0 - <j> + (j3 - a) A) + e((f> + A)) . 
Adding w 2 e (9 + X + (j3 — a) A) to both sides gives 

e (0 + A + 03 - a) A) = -1 + w (e (0 - <j> + (/? - a) A) + e (0 + A)) . 
Subtracting u> e (0 — + (f3 — a)X) yields 

e (0 + A + (jS - a) A) - w e (0 - </> + (0 - a) A) = -1 + iw e (0 + A) . 
Factoring the left hand side 

(e (A) - to e (-0)) e (0 + (/? - a) X) = -1 + w e {<j> + A) 

or 

(e + A) - w) e (0 - <f> + (0 - a) A) = -1 + w e {4> + A) . (3.19) 
Rearranging (3.19) gives (3.13). ■ 

The following lemma established some simple facts about the right hand side of 
(3.13). 

Lemma 3.12. Fix < w < 1. Consider the Mobius transformation 

11)7 — 1 

Mz := , zeC. (3.20) 

z — w 

Then M = M w is a bijection of the unit circle T onto itself and there is a decreasing 
continuous bijection g of the real line M onto itself such that g(0) = —1/2 and 
Me(t) = e (g (t)) for all t G E. 

Proof. Recall any Mobius transformation is a continuous bijection of the Riemann 
sphere onto itself. Note M is its own inverse. Taking the modulus we see 

|iue(i)-l| = |e(t)-u;|, t G M. 

Hence M maps T into itself and consequently, M is a bijection of T onto itself. 

In particular, ^(t) := Me(t), t G R has minimal period equal to one and maps 
any interval of non-zero integer length onto T. Since M has positive coefficients M 
maps the extended real line onto itself. Since Mi is in the first quadrant M maps 
upper half-plane onto itself, in particular, 7 ([0, 1/2]) is the part of T in the upper 
half-plane. Note 7 (0) = Ml = -1 = e (-1/2) . 

By the path lifting theorem there is a unique continuous function / : [0, 1] — >• M 
such that /(0) = -1/2 and e (/(f)) = j(t) for t G [0, 1]. Since e (/(1/2)) =7(1/2) = 
M(-l) = 1 and e (/ ([0, 1/2])) = 7 ([0, 1/2]) is the part of T in the upper half-plane, 
it follows that /(1/2) = —1. Consequently, / is a continuous bijection of [0, 1] onto 
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The desired function is determined by 



g(t + n) := f(t) - n 



for < t < 1 and integers n. 



Remark 3.13. Since ff(l/4) is in the first quadrant and ff(l/2) = 1 then e (g [0, 1/4]) 
is the part of T in the fourth quadrant and part of the first quadrant and e (g [1/4,1/2]) 
is the remaining part of T in the first quadrant. This "explains" the general shape 
of the curve in our pictures below. 

We are now ready to finish our analysis of the spectrum of Pb- 

Theorem 3.14 (0 < w < 1). Let the two intervals Ii and I2 be as described, with 
a and (3 fixed and denoting the endpoints of I2, i.e., 1 < a < (3. Pick parameters 
4>, if), and w, < w < 1, describing one of the self adjoint extensions Pb- Then the 
spectrum of Pb is pure point, with the eigenvalue list equal to the sequence {A n } of 
real numbers solving equation (3.13). Every number \ n in the spectrum of Pb is 
isolated, and, for n fixed, the corresponding eigenspace is spanned by the associated 
eigenfunction from (3.15). 

Proof. By Lemma 3.11 the eigenvalues of Pb are the real solution to 

e (0 - <t> + {/3 - a) A) = M e (cp + A) 

where M is the Mobius transformation (3.20). Hence, if g is the function from the 
conclusion of Lemma 3.12, then the eigenvalues are the solution to 

e (6 - <j) + {P - a) A - g(<f> + A)) = 1. 

Hence the eigenvalues are the A for which 

6-<t,+ (P-u)\-g{<t) + X) SZ. (3.21) 

Let h : M — >• R be determined by 

h(t) :=9-<t>+(f3-a)t-g(cf) + t). (3.22) 

If t is an integer and n is an integer, then g(t + n) = g(t) — n, consequently, 

h ([i, t + n)) = [h(t),h(t) + (1 + [ft - a)) n) . (3.23) 

is an interval of length n(l + /3 — a). In particular, it contains at most 1 + 
n (1 + (3 — a) integers, hence n points from the spectrum of Pb- B 



Figure 3.1 shows three plots of the left hand side of (3.21) as a function of t, the 
plots illustrate that the curve is less linear for larger values of w. 
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Figure 3.1. y = 9 - <f> + (/? - a)t - g(<p + t) for 6,i/>,<f> = 0, 
p - a = 1, and w = 0.1, 0.5, 0.9. 

The spectrum of Pb is continuous in B in the sense that 

Corollary 3.15. For integers n, let X n denote the solution to 8 — <f> + (/3 — a)t — 
g(4> + t) — n. The the spectrum of Pb is the set {X n | n G Z} , A„ < A„+i, and for 
fixed n, (w,(f),^,8,a, /3) — > X n is continuous. The dependence of the eigenvalues on 
the parameter w in B is harmonic. Specifically in the formula (3.24) the circle is 
given by the parameter t, and when viewed as a function on the circle, we therefore 
get an extension to the inside disk, represented in polar coordinates (w,u) where w 
is radius and u is angle; and with continuity up to w = 1. 

Proof. Let h be as in (3.22). Then A„ is the first coordinate of the point of inter- 
section of the curve y = h(t) and the horizontal line y = n. Clearly the curve is 
continuous in (f>, ip, 9, a, and j3. By construction, g is the argument of a continuous 
logarithm. Consequently, 



where 7(i) — Me(t) as in Lemma 3.12. To get the harmonic extension, note that 
the Poisson kernel in (3.24) offers a harmonic extension from the circle to the inside 
disk. We use the parameter t for the circle, and the extension is to the inside disk 
is in polar coordinates (w, u) with w as radius and u angle. So the statement about 
continuity up to w = 1 is just continuity up to the boundary circle for the harmonic 
extension. The latter is a known property of Poisson kernel. ■ 

Remark 3.16. The spectrum of Pb is uniform in the sense that for any integer n, 
we have 

h ([An, A„ + l)) = [n,n + l + P-a). 

Where A n is as in Corollary 3.15. Consequently, the interval [A„, A n + 1) contains 
exactly [1 + f3 — a\ points from the spectrum. Here, [1 + /3 — a\ is an integer such 
that f3-a < [1 + (3 - a\ < I + [3 - a. 

Points in the spectrum are separated in the sense that 

Corollary 3.17. There exists 5 > 0, so that if X ^ A' are eigenvalues of Pb, then 




(3.24) 



| A — A' 



\>5. 
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Proof. The function t — > (e t | 1) = J Q e t + f a e t is continuous and = 1 at t = 0. 
Hence there is a 5 > 0, such that \t\ < 5 implies \(e t | 1)| > \- Since A — A' is a root 
of t -> (e t | 1). we conclude |A — A'| > <5. ■ 

Let A„ be the enumeration of the spectrum of P from Corollary 3.15. The 
following figures shows plots of A ra as a function of w for several values of n. Here 
both the selfadjoint operator Pb, and the numbers A„ in its eigenvalue list, depend 
on the choice of the matrix B, but for brevity, in some formulas, we have suppressed 
the B-dependence. Now A„ = X n (B), but in Figure 3.2-3.4 we zoom in on the 
dependence of the w parameter from B. Figure 3.2-3.4 show examples where S w —> 
as w — > 1 and where S w > S > for all w. Here 6 W = S is determined in Corollary 
3.17 above. 



Figure 3.2. ft = (0, 1) U (2, 3), and 6, ip, 4> = 0; A w=0 = \ + \l, 
and A„, = i = Z. 



Figure 3.3. ft = (0, 1) U (2, 4), and 6, ip, 4> = 0; A w=0 = \ + |Z, 
and A„, = i = ^Z. 



FIGURE 3.4. il = (0,1) U (2,3), ^0 = 0, and 9 = 1/2; A w = |Z, 
for all < u> < 1. 
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Corollary 3.18. Suppose j3 — a is rational and p, q > are integers such that 
P — a = p/q. Then there is a finite set L such that the spectrum of Pg is L + qZ. 

Proof. Suppose p > 0. Let A„ be as in Corollary 3.15. In particular, h (Ao) = and 

h (A + q) = h (A ) + (1 + /3 - a) q = q + p. 

Hence, setting t = Ao and n = q in (3.23) gives 

h([\ ,X +q)) = [0,q + p) 

which contains q + p integers. Let L :— {Ao, Ai, • • ■ , A g+p _i} be the corresponding 
eigenvalues of Pg. By construction L C [Ao,Ao + <7) and the spectrum of Pg is 
L + qL. M 

3.3. Irrational f3 — a. In the remainder of this section we discuss the spectrum of 
Pg when ft — a is irrational. 

Corollary 3.19 (Asymptotics). Let A„ be the enumeration of the spectrum of P 
from Corollary 3.15. If f3 — a is irrational, then 

A + — 1 < A fe < A + ^ X + 1 

1 + p — a 1 + p — a 

for all k £ 1. 

Proof. For any real number x, let \x\ be the real number satisfying < \x\ < 1 
and x — \x\ £ 7L. Well, for any integer q > 0, the set h ([Ao, Ao + <?)) is the in- 
terval [0, (1 + /? — a) q). By irrationality of (3 — a, the interval [Ao,Ao + (/) con- 
tains [(1 + /3 — a) q] points from the spectrum. Similarly, [A , Ao + q + 1) contains 
[(1 + j3 — a) (q + 1)] points from the spectrum. Consequently, 

Ao + q < A[(i +( 3_ Q )q] + i < Ap(i +( g_ Q ) g ] +2 < • ■ • < Ap(i +( g_ Q )( g+ i)] < Ao + q + 1. 

Restating we see 

A +<z< A fc < Ao + 9+1, (3.25) 
when [(1 + P - a) q] + 1 < k < [(1 + /3 - a) (q + I)] . Now 

[(1 + - a) q] + 1 < k < \(l + (3 - a) (q + 1)] 

implies 

(l + p-a)q + l<k<(l + p-a)(q + l) + l. 
Solving for q gives 

fc-1 

q< t— 5 <«+!■ 

1 + p — a 

The desired conclusions now follow from (3.25). 

We leave the discussion of the case q < for the reader. ■ 

In the general case (3.14) we study the intervals in the A axis between neighboring 
asymptotes: a sequence of equal-length intervals. The difference p — a: rational vs 
irrational. We now study the points on the axis corresponding to the intersections 
in the above examples. By Lemma 3.11 these points constitute the spectrum A. 

Remark 3.20. If p — a is irrational, the spectrum is purely aperiodic, see details 
below in Example 3.21. 
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The examples below serve to illustrate a number of spectral theoretic issues 
accounted for in our theorems. 

They are motivated by several considerations: 

(i) The examples, both in this section, and in Section 5, show that the possibil- 
ities, outlined in the abstract, may indeed arise. 

(ii) A vector bundle over U(2): This is illustrated below. Recall, for each B in 
U(2), the eigenspaces for Pg form a line bundle over the group U(2) as base-space, 
with degeneracies at places when the multiplicity jumps from 1 to 2; 

(iii) The examples further serve to make concrete the subtle interplay between 
geometry and spectrum; 

(iv) They illustrate some additional issues dealing with the way a continuum 
of parameters has symmetry breaking, and thus produce cases with properties of 
special interest; for example illustrating how the spectral pairs (SP) form a small 
subfamily; and 

(v) The various families of two-interval cases serve as models for scattering theory 
in physics (see Section 2.12.1), so they are models for systems with a much richer 
structure, models for potential barriers, for quantum jumps, ... . By model, we 
mean, up to unitary equivalence, component-by component. 

Example 3.21. Choose w = 9 = -1/4, <f> = -1/8, a = 3, (3 = 3 + y/2. 

Let J s = [s + 1/8, s + 1 + 1/8), s £ Z, i.e., J s is the s th interval between two 
neighboring branch cuts. Let L s be the eigenvalues in J Sl then A = {J seZ L Sl and 
the set U sS z (L s — s) is dense in Jo. See Figures 3.5-3.7. 

This conclusion holds whenever I2 is chosen with irrational length {3 — at. The 
reason is that this number then induces an irrational rotation which is known to 
have dense orbits; a basic example in symbolic dynamics. 

Figures 3.5-3.7 are based on (3.13), the argument function from the proof of 
Lemma 3.12, reduction modulo one gives the branch curves. More precisely, the 
figures are obtained in a number of steps. Step 1. Fix B as in (2.10), and let 
Pb be the corresponding selfadjoint operator. Consider the curves for argument 
function from the right hand side of (3.11) and the lines with slope j3 — a from 
the left hand side of (3.11). We are assuming here that this slope is irrational. 
Step 2. Identify the asymptotes, and fix an interval between two branch cuts; for 
example the interval Jq on the axis closest to 0. Note that all the intervals between 
neighboring branch cuts have the same fixed unit-length. They extend both to the 
left and to the right of Jo- Step 3. Since the spectrum of Pb is discrete and 
infinite, it intersects all these intervals between branch cuts. Now, translate all of 
these finite intersections down to Jo. Conclusion. Since the line-slope is irrational, 
the union of all these sets (inside Jo) is dense in Jo- 

Remark 3.22. When the difference /3— a is irrational, then the distribution of points 
in A changes as you move around on the axis between the cut-intervals. See Remark 
3.20, and Example 3.21 above for a detailed discussion. 

Corollary 3.23. Let X n be the enumeration of the spectrum of P from Corollary 
3.15. For a real number x, let [x] be the real number satisfying < [x] < 1. If (3 — a 
is irrational, then [X n ] ^ [A m ] for all n ^ m and the set {[A„] | n € Z} is dense in 
the interval [0, 1]. 
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Lo UI =0 (L S - s) 



Figure 3.5. Dense orbits 




U 2 s=0 (L s -s) U 3 s=0 (L s -s) 



Figure 3.6. Dense orbits 




U 7 s=0 (L s -s) \Jll a (L s -s) 
Figure 3.7. Dense orbits: The spectrum as a random-number generator. 

Proof. A„ is a solution, for t, to 

0-(j>+(P-a)t = n + g(<f> + t) (3.26) 

iff 4> + A n is a solution, for t, to 

6-(l + P-a)<j>+(J}-a)t = n + g(t). 

Consequently, replacing by corresponds to translating the spectrum and replac- 
ing 6 by — (1 + j3 — a)4>. We will therefore assume (j> = in (3.26). 

Suppose [A m ] = [A„] , Then A m = A n + k for some integer k. Hence M e (A m ) = 
Me (A„) and therefore Lemma 3.11 implies e (9 + (f3 — a) A m ) = e (8 + (0 — a) A„) . 
Consequently, 

9 + (/3-a) (A„ + k) = 6 + {fi - a) A„ + j 

for some integer j. It follows that f3 — a= I . Contradicting that /3 — a is irrational. 

We proceed the proof of the density claim by an alternative description of the 
set {[A„] | n e Z}. For real x let — | < {x} < | be chosen such that x — {x} is an 



MOMENTUM OPERATORS IN TWO INTERVALS 



19 



integer. Let / be the path lifting function from the proof of Lemma 3.12. Replacing 
t in (3.26) by t + m where < t < 1, we get 

{8 + (j3-a)(t + m)} = l + f(t) 

since {n + g(t + m)} = {n + /(f) - m} = l+/(t). Let l m (t) := {6 + (/? - a)(f + m)} . 
We will say £ m is a line. Then, the set {[A„] | n € Z} is the set of solutions to the 
equations 

*m(*) = 1 + /(*), 

in the interval [0, 1). 

To show that {[A„] | n G Z} is dense in the interval [0, 1] we must show that any 
subinterval [a, b] C [0, 1] contain at least one point of the form [A„] . Consider the 
rectangle R = [a,b] x [1 + /(&), 1 + f(a)]. Note the graph of 1 + /(f), a < t < b 
is a monotonely decreasing function passing through the upper left hand corner 
of R and the lower right hand corner of R. Since {[m(f3 — a)] m e Z} is dense 
in [0,1], there are lines £ m passing through a dense set of point of the bottom 
[a, b] x [1 + f(b), 1 + /(&)] of R. Since the lines £ m all have fixed positive slope j3 — a 
we can ensure that they pass through the right hand edge [b, b] x [1 + fib), 1 + /(a)] 
of R. These lines must intersect the graph of 1 + /(f), the first coordinates of these 
intersection points are the points in { [A„] | n € Z} . ■ 

Proposition 3.24. Let Ji,i2 &e the two intervals as before, and Pb one of the 

self adjoint extensions. Let A be the spectrum of Pb- Then A = F + Z, where F is 
some finite subset in M., if and only if ft — a is rational. 

4. Unitary One-Parameter Groups 

Consider the Hilbert space ,34? — L 2 (il), Q, = Ji U I2, and with the two intervals 
Ji and I2 as described. As before let 

\ VI - w 2 e(V>) we(6-(j>) ) 1 ' ' 

be a unitary matrix with the parameters w, ip, 9, ip chosen. 

Consider the Hermitian operator L — ^ ^ on its minimal domain ^0 corre- 
sponding to function / € such that /' e 3tf and /(0) = /(l) = f(a) = f{(5) = 0. 
Then the selfadjoint extensions Pb are prescribed by the choice in (4.1). Specifi- 
cally, 

dom(P B ) = {/ G dom(P) \f R = Bf L ] (4.2) 
where P is the maximal operator and where 

We get a unitary one-parameter group Us{t) '■ ffl — > Jrf?, ( 6 K such that for 
/ € dom(P A/ ) 

Km~(U B {t)f-f)=iP B f. (4.4) 
The following is a consequence of the results in Sections 2-3. 
Corollary 4.1. 

(1) Let f be some wave-function localized in Iy, then if x and x — t are in l\, 
then 

(U B (t)f) (x) = f(x-t). 
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(2) As the support ofU B (t)f hits x = 1 (the right-hand side boundary point), 
then it transfers to the two parts (left-hand side point in Ii), and to a 
(left-hand side part in I2) with probabilities w 2 and 1 — w 2 , respectively. 

(3) In the transfer from x = 1 to x — 0, the phase is shifted by e{(f); and from 
x = 1 to x = a, it is phase-shifted by e(ip). 

(4) The boundary conditions (4-2) are preserved by U B (t) for all t € R; i.e., 
we have 

(U B (t)f)(l) \_ f (U B (t)f)(0) \ 
(U B (t)f)(P) J *\{U B {t)f){a) ) W 

for all f £ doTa(P B ); see also eq. (4-3). 

Proof. It follows from (4.4) that dom(P B ) is invariant under U B {t) for all t G M, 
and (4.5) follows from this and Corollary 2.3. ■ 







FIGURE 4.1. Quantum Binary Model 



5. Spectral Pairs 

In the study of the spectral pair-problem for subsets in M. d ([Pug74, J0r81, JP98, 
JP99]) one considers subsets fl in K d of finite positive measure. If the Hilbert space 
L 2 (il) has a Fourier basis, we say that is a spectral set. This means that there is 
a discrete subset A in M. d such that E(A) :— {e\ | A G A} is an orthogonal basis in 
L 2 (il). We then say that (O, A) is a spectral pair, and the set A is called a Fourier 
spectrum. The spectrum is never unique. We say £1 is a tile, or that Q tiles by 
translations, if there is a set A such that J2 a eA Xn(t — a) = 1 for a.e. t € R d . The 
set A is called a tiling set for Q, 
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An old question from [Fug74] asks if SI is spectral if and only if it tiles R d by 
translations. The question is known to have negative answers for d > 3. But it is 
open for d = 2, and 1. Another reason for our focus here on d = 1 and Q = union 
of two intervals, is that, in this case, we prove in Corollary 5.13 below that the 
two properties are then equivalent. In other words, Fuglede is affirmative in ID for 
these sets. 

While the equivalence of the two notions, "tile by translations" and "spectral 
set" is open for general measurable subsets SI in R, of finite positive measure, we 
show below that, in case when S7 is a union of two intervals, then the equivalence 
is affirmative. For a discussion of the question for d > 3, see e.g., [KM06, IKT03, 
Lab02, IP98, Ped96, LW97]. The first result proving non-equivalence was for d = 5, 
the paper [Tao04]. Others followed shortly after, proving non-equivalence for d = 4, 
and d = 3. 

Let 1 < a < p and suppose Pb is the selfadjoint restriction of the maximal 
operator P corresponding to the matrix B — B(w, (f>,ip,9) in (2.10). Let As be the 
spectrum of Pg. By (3.1) any eigenfunction of Pb is of the form 



for some complex numbers a\ and b\. We say Pb is spectral if a \ = b\ for all 
A £ A B . 

Lemma 5.1. If Pb is spectral, then ([0, 1] U [a, f3] , Ab) is a spectral pair. Con- 
versely, if ([0, 1] U [a, 0\ , A) is a spectral pair, then there is a spectral Pb such that 
A A B . 

Proof. Since the eigenfunctions of Pb form an orthogonal basis the first claim is 
obvious. 

Conversely, suppose ([0, 1] U [a, 0], A) is a spectral pair. Then the closure of 



determines a selfadjoint restriction P of P. The spectrum of this selfadjoint restric- 
tion is A and the eigenfunctions are (the multiples of) the ca's with A <E A. Since any 
selfadjoint restriction of P is of the form Pg for some B the proof is complete. I 

When ([0, 1] U [a,/3],A) is a spectral pair eq. (3.3) can be solved to find a 
corresponding matrix B. 

Theorem 5.2. Let Ab be the spectrum of the selfadjoint restriction Pb corre- 
sponding to the matrix B = B{w,<j),^p,9) in (2.10). If w = 1, then Pb is not 
spectral and consequently ([0, 1] U [a,/?], As) is not a spectral pair for any choice of 
the remaining parameters a, f3, 9, if), tfi, w. 

Proof. By Theorem 3.3 every eigenfunction has a\ = or b\ = 0. Hence this is a 
direct consequence of Lemma 5.1. I 

5.1. Two Theorems. With a choice of two intervals, specified by the endpoints 
a and j3 for the second interval, we proceed to characterize the four parameters in 
the boundary matrix B from (2.10) in Corollary 2.3 which yield a spectral pair. In 
other words, we determine when there is an orthogonal Fourier basis in I? of the 
two intervals. Our results are Theorems 5.3 and 5.8 below. Together, they exhaust 
the possibilities for spectral pairs formed from the union of two intervals. 
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In particular, we show that all the two-interval spectral pairs fall in two classes, 
(i) lattice tilings (Theorems 5.3), and (ii) non-lattice tilings (Theorems 5.8). In both 
cases, we find the particular instances of spectrum(Pe) occurring as the second part 
in a spectral pair. Since we are in one dimension, (i) represents the cases when the 
two intervals tile R by a fixed multiple of the integers Z. In this Fourier 
spectrum may be taken to be the dual lattice in E. We show that, for (ii), the 
possible Fourier spectra arising as spectrum(Pe) must be the union of two lattice 
co-sets. 

Theorem 5.3. Let Kb be the spectrum of the selfadjoint restriction Pb corre- 
sponding to the matrix B = B(w,^),tf>,0) in (2.10). If w = 0, then Pb is spectral 
iff 

eN (5.1) 



1 + p-a 



and the parameters <j>, 8 satisfy 



(0-!)(l-a) 



In the affirmative case the spectrum is As = \+p- a + 
Proof. By Theorem 3.9 and Lemma 5.1 Pg is spectral iff 



for all ri e Z. That is iff 



l-e+n\ 



(6- i) (l-a) + /3n 
1 + p — a 



for all n £ Z. Considering n = and n = 1 it follows that Pb is spectral iff the 
conditions (5.2) are satisfied. 

The formula for the spectrum is from Theorem 3.9. I 

Corollary 5.4. If (5.1) and (5.2) then ^[0,1] U [«,/?], 1+ ^_ a Z^ is a spectral pair 
and [0, 1] U [a, j3] is a tile with tiling set (1 + (3 — a)Z. 

Proof. Iik= 1+ l_ a , then 

[a - (1 + p - a)(k - l),/3 - (1 + p - a)(k -l)] = [l,l+0- a], 

consequently, [0,1] U [a, j3] is a tile with tiling set (1 + j3 — a)Z, when k is an 
integer. ■ 

Remark 5.5. In Theorem 5.3 the condition (5.1) restricts the geometry of the in- 
tervals. Specific examples of sets satisfying this condition includes [0,1] U [|,3] 
and [0, 1] U [7ir + 1, 14ir] . The other condition in (5.2) is a restriction on the set of 
selfadjoint restrictions, that is on the values of 9 and ip, that determine B and the 
spectrum. 

For every j3 — a £ Z + , we may construct tiling spectral pairs. In fact, fix 
p = (3 — a, and set p — 2(1 +p), a = 2 + p. Then 

O = [0,l]U[2+p,2p + 2] 

tiles with (1 + p)Z. 
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Example 5.6. Suppose j3 — a = 2. The condition 1+ ^_ a = § € Z, see (5.2), 
implies /3 e 3Z. For ,3 = 6, the set [0, 1] U [4, 6] tiles K with 3Z; see Figure 5.1. By 
contrast, [0, 1] U [2,4] does not tile, and note that (5.2) does not hold. 



Figure 5.1. SI U (ft + 3) U (ft - 3) U (ft - 6), where SI = [0,1] U [4,6]. 

Example 5.7. Suppose (3 - a = 3. Then 1+ ^_ a = f € Z implies /3 e 4Z. For 
(3 = 8, the set [0, 1] U [5, 8] tiles with 4Z, by Theorem 5.3. 



FIGURE 5.2. QU (ft + 4) U (ft -4), where ft = [0,1] U [5,8]. 

Theorem 5.8. Let Ab be the spectrum of the selfadjoint restriction Pb correspond- 
ing to the matrix B = B(w,<j),ip,ff). If < w < 1, then Pb is spectral iff 

1 < a is an integer, /3 = a + 1, (5-3) 

and the parameters 9,(j),ip, and w satisfies 



k e Z }. (5.4) 



1 f / 1 + 2k 

9 — 2(f) e Z, ip + (a — 1) <f> E -Z, and w € < cos 27r 

2 [ \ 4a 

In tfte affirmative case Ab = | — ± cos 27r ^ ul ^ | + Z, hence Ab = { —0 ± ^ 2fc } + Z, 
for some integer k. 

Proof. We saw earlier that A is a eigenvalue for Pb iff 

e(0-< t > + ( f 3-a)A) = We V + V- 1 (5.5) 
e (0 + A) — w 

and if A is an eigenvalue, then a\ = b\ = 1 iff 



VI - e + A - aA) = 1 - w e (0 - + (/? - a) A) . (5.6) 

Taking the square of the modulus to both sides of (5.6) we get 

= 2w (w - c (6 - + (p - a) A)) . 

Where c(t) := cos(2-7rf). Since w ^ 0, we conclude 

w = c(9- <j>+ (f3- a) A) . (5.7) 

Replacing w in (5.5) by c (9 — cf> + (f3 — a) A) and simplifying using s(t) := sin(27ri) 
and e(t) = c(t) + is(t) we arrive at 

i (c (<j> + A) - c (9 - <$> + ((3 - a) A)) s (9 - $ + (f3 - a) A) 

= -1 + c 2 (6» - + (0 - a) A) + s (4> + A) s (9 - cf> + ([3 - a) A) . 
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Consequently, 

(c (0 + A) - c (0 - (}> + (/? - a) A)) s (0 - <f> + {0 - a) A) = (5.8) 

c 2 (0 - <t> + (/3 - a) A) + s (4> + A) s (9 - <j) + (0 - a) A) = 1. (5.9) 

Since < w < 1, if follows from (5.7) that s (0 - + (/? - a) A) 7^ 0, hence (5.8) 
implies 

c {(j) + A) = c (9 - (f> + ((3 - a) A) = 10. 
It follows from (5.7) and (5.9) that 

s ((f> + A) = s (9 - cf> + ((3 - a) A) = ±^1 - w 2 . 

Consequently, 

e (0 - + (13 - a) A) = e (0 + A) = w ± i\A - w 2 . (5.10) 



is one of the two fixed points w ± iyl — w 2 of the Mobius transformation M = M w 
from Lemma 3.12. Note (5.10) implies (5.5). Hence, if A satisfies (5.10), then A is 
in K B - Plugging (5.10) into (5.6) we get 

e ("0 + A — a\) — \/l — w 2 ± iw — ±i (w =p i\J 1 — w 2 ^j . (5-H) 

We have seen that A satisfies (5.5) and (5.6) iff A satisfies (5.10) and (5.11). 
Let A := cos -1 (u;)/27r. Then 



e (i-M — w ± i\J 1 — w 2 

by (5.10). And c (4> + A) = w, implies A g |— <fi ± Aj + Z, hence 

A s C |-</»± A} +Z. (5.12) 
By [Lan67] As has density 1 + j3 — a > 1. Hence one of 

K B n ({-c/> + A j + z) and A B n ({-0 - A j + z) 

has density > 1/2. Assume it is K B n + ^} + Z J • 

Let be the integers fc such that A^ := —0 + A + k £ As and similarly, let Z_ 
be the integers k such that A^T := — — A + k £ A#. Plugging A = A^ = —0 + A + k 
into (5.11) gives 

e (V - (</» + A) (1 - a) + (1 - a) kj = ie = e Q - A^ . 

for fc G Z + . Rewriting we get 

e ((1 -a)A0=eQ + aA^j e (0 (1 - a) - ip) 

for k £ Z + . Since the right hand side is independent of k £ Z + and Z + has density 
> 1/2, 1 — a is an integer. Similarly, it follows from (5.10) that (3 — a is an 

integer. The set |— </>±a| + Z has density equal to two, hence the subset Kb 

has density at most two, so by [Lan67] 1 + /3 — a < 2. It follows that (3 = 1 + a 

and As = |— </> ± a| + Z. In particular, Z+ = Z_ = Z, a is an integer > 1 and 

/3 = 1 + a. 
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Using /3 — a = 1 and A 



— <f> ± A + k for k € Z we can write (5.10) as 




thus — 20 is an integer. Similarly, (5.11) becomes 



e (V> + (a - 1) 0) 



= e 



( 




Consequently, 



±- ± aA G -0 + 



1)0 + Z. 



Adding these equations shows 2 + (a — 1) (j)) € Z. Subtracting them gives 2a A € 



Corollary 5.9. If (5.3) and (54), then (JO, 1] U [a,/3], |± cos 27r w j + Zj is a spec- 
£raZ pair and [0, 1] U [a, /3] is a tile with tiling set {0, 1, 2, . . . , a — 1} + 2aZ. 

Remark 5.10. In Theorem 5.8 the condition (5.3) restricts the geometry of the in- 
tervals. The set [0, 1] U [2, 3] satisfies (5.3) but not (5.1). The set [0, 1] U [7?r + 1, 14tt] 
satisfies (5.1) but not (5.3). Finally, [0, 1] U [3,4] satisfies both (5.1) and (5.3). The 
other condition in Theorem 5.8 (5.4) is a restriction on the set of selfadjoint restric- 
tions, that is on the values of w, 8, </> and ip, that determine B and the spectrum. 

Remark 5.11. (i) Let il :— [0, 1] U [a, /?]. Suppose (5.1) or (5.3) is satisfied, then 
is a spectral set. Any choice of a unitary matrix gives us a spectrum A#. Most of 
these spectra will not satisfy (5.2) or (5.4). Hence for most B the pair (f2,A#) is 
not a spectral pair, even when VL is a spectral set. 

(ii) Suppose (f2, A) is a spectral pair. By Lemma 5.1 there is a unitary matrix 
B, such that Pb is spectral and As = A. Construct a new unitary matrix Bq by 
replacing ip by a -0o, (keeping w,<p, and 9 fixed) such that the appropriate one of 
(5.2) or (5.4) fail. Then Pg is not spectral, yet, since the spectrum is independent 
of ip we have Ab = As = A. 

(iii) If (f2, A) is a spectral pair we can choose selfadjoint restrictions P\ and P2 
of the maximal operator P such that A = Ap 1 = Ap 2 and the functions e\, A € A 
are eigenfunctions for P± but not for Pi. 

Corollary 5.12. The set [0, l]U[a, 0\ is a spectral set iff either 1+ p_ a is an integer, 
or a > 1 is an integer and (3 = a + 1. 

Corollary 5.13. // the set [0, 1] U [ot,0\ is a spectral set, then it tiles the real line 
by translations. 

Corollary 5.14. If [0, 1] U [a, j3] is a spectral set and f3 — a > 1, then [0, 1] U [a, f3] 
is a lattice tile. 



2 + Z. But this is equivalent to 




This completes the proof. 
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5.2. Applications. While our results for spectral pairs built with two intervals 
offer the list of possibilities in terms of the two parameters a, /?, and the boundary 
matrix (2.10), it is helpful to visualize some of the possibilities, and discuss geomet- 
ric implications. For example, even though the sets spectrum(Pe) is determined 
in all cases, its properties are reflected very concretely by curve intersections, see 
Figures 5.6 through 5.12 below. And our brief discussions of the examples below 
serve to highlight main points. For example, one arrives at a visual for the fact that 
the difference j3 — a represents a slope. Also when it is irrational the restriction 
to spectrum(Ps) of natural homomorphism from K to the quotient E/Z (the circle 
group) has dense range. See Figures 5.4 and 3.5. 

In the discussion below, we further illustrate how the possibility for spectral type 
for our 6-parameter family of selfadjoint operators depend the roots of modulus 1 
in certain characteristic polynomials; see Example 5.18, Figure 5.8; Example 5.19, 
Figure 5.9; and Example 5.21, Figure 5.11. In detail; depending on the relative 
position of two intervals with integer endpoints, we get a complex polynomial P(z). 
The question of how many roots are on the circle \z\ = 1 decide a number of 
properties of the spectrum of an associated operator Pb- 

Example 5.15. Let Q = [0, 1] U [a, /?], and fix the matrix B = B(w, cj>, t/j, 8). The 
spectrum of Pb is determined by 

e (<p + a) — w 

See Lemma 3.11 and 3.12 for details. That is, A E As iff it is at the intersection of 
the linear function 9 — <f)+((3 — a) A (mod 1), and the monotone decreasing function 
g (t) (mod 1). Density of As is reflected in the slope (f3 — a). 

Choose w — 9 = —1/4, </> = —1/8, a = 3, and let (3 = 4,5,6,7. Since 
j3 — a is rational, the spectrum is periodic; sec Figure 5.3. On the other hand, 
let P = 3 + {1, 2, 3, 4}y/2, so that f3 — a is irrational, and the spectrum is purely 
aperiodic; see Figure 5.4. 



5.3. Examples. What follows are examples 5.16 - 5.22; and a summary of selected 
conclusions from these examples. 

In Examples 5.16 and 5.17, we compute the boundary matrix B for the case 
when a = 2, = 3. The union of the corresponding two intervals offers a spectral 
pair which is not a lattice tile. 

Continuing with cases when the second interval has the form [a, ft] with integer 
endpoints, we follow with a family of examples /3 — a = 1; i.e., same length. They 
again are spectral pairs, but when a is odd, then they are also in the more restrictive 
class, the lattice tiles. Indeed, for a odd, it turns out that there is then both a 
lattice spectrum, and a non-lattice spectrum. Sec Example 5.20 which also shows 
that different choices of boundary matrix B, while leading to the same spectrum, 
yet can result in different systems of eigenfunctions. In Figures 5.7 - 5.10, we 
illustrate determination of spectra from curve crossings. Examples 5.21 - 5.22 
illustrate computation of spectrum(Pe) for sets that are unions of pairs of intervals, 
different length, and producing non-spectral pairs, i.e., not spectral pairs for any 
choice of A. 
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Example 5.16. We know that [0, 1] U [2, 3] is a spectral set with spectrum Z U 
(| + Z) . In this case B = ( ™ ^ ~\ must satisfy the equations 



::)(:)=(! 

w x \ I i \ _ ( 1 

y * A -* A V -i 



A integer 

, A integer plus 1/4 



Hence 



The boundary matrix B in (5.13) has parameters w — v2/2, <^> = — |, ^ = g, 
and 8 = —4. By (3.13), A is in the spectrum if and only if 

1. n/2 e(A-|)-y / 2 

Since both sides of the above equation have modulus one, comparing arguments 
shows that 

i.e., A = {0, |} + Z = ZU(| +Z) confirming the claim at the beginning. We revisit 
this example several times below. 

In the following examples we apply Corollary 5.12 in order to decide when L 2 (JiU 
I2) has a Fourier basis, i.e., when we get spectral pairs. In the examples below this 
amounts to counting the number of root on the unit circle of a certain polynomial, 
see also [Ped96]. 

We illustrate the results with two classes of examples, see Figures 5.5, 5.6. Fix 
the parameters by w = 6 = — 1/4, <j> = — 1/8, i\) = 1/8, and consider the two 
intervals I lt J 2 given by [0, 1] U [2, 2 + k] in Class A, and [0, 1] U [3, 3 + k] in Class 
B, where fc = 1, 2, 3 . . .. Note that 

(1) The spectrum A depends on j3 — a but not the location of the interval [a, f3]; 

(2) The coefficient a(A) depends on all 6 parameters. 






1 2 
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Figure 5.5. Class A: [0, 1] U [2, 2 + fc] 
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Figure 5.6. 


Class B: [0,1] U [3,3 + fc] 



Example 5.17. Spectral pair: [0,1] U [2,3] , w = 6 = -1/4, <p = 

V = 1/8. 
Spectrum 

A = {0, 1/4} + Z = Z U (1/4 + Z) 

Coefficients 

a A = 1,VA € A 




p(z) = (z- l)(l + z 2 ) A 

Figure 5.7. 



Example 5.18. Not spectral pair: [0, 1] U [2,4] , w = 4g, Q = -1/4, (f> 

ip = 1/8. 

Spectrum 

Aw {0,0.1825, 0.5675} + Z 

Coefficients 



ao 

00.1825 
°0.5675 



1 

1.61716-0.455719? 
-0.367157 + 0.205719i 
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p(z) = (z-l)(l + z 2 (l + z)) A 



Figure 5.8. 

Example 5.19. Not spectral pair: [0, 1] U [2,5] , w = = -1/4, (f> = -1/8, 

ip = 1/8. 
Spectrum 

A « {0, 0.1550, 0.3913, 0.7037} + Z 

Coefficients 

a = 1 

ao.1550 « 2.07306 - 0.464978z 

ao.3913 « 0.253715 + 0.489548z 

ao.7037 « -0.326779 - 0.27457i 




p(z) = (z- l)(l + z 2 (l + z + z 2 )) A 

Figure 5.9. 

Example 5.20. Spectral pair: [0,1] U [3,4] , w = -2=, = -1/4, cf> = -1/8, 

V = 1/8. 
Spectrum 

A = {0, 1/4} + Z = Z U (1/4 + Z) 

But the coefficients are 

«o = 1, ai/4 = i 

so the parameters above lead to non-spectral eigenfunctions. Illustrating that dif- 
ferent parameters B = B{6, ip, cf>, w) can lead to the same spectrum yet to different 
eigenfunctions. 
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p(z) = (z- l)(l + z 3 ) A 



Figure 5.10. 

Example 5.21. Not spectral pair: [0, 1] U [3, 5] , w = 6 = -1/4, = -1/8, 

$ = 1/8. 

Spectrum 

A « {0,0.1825,0.5675} +Z 

Coefficients 

a = 1 
ao.i825 « 1.08072 + 1.28644* 
ao.5675 « 0.419281 - 0.0364378i 




p{z) = (z-l)(l + z 3 (l + z)) A 



Figure 5.11. 

Example 5.22. Not spectral pair: [0, 1] U [3, 6] , w = 6 = -1/4, = -1/8, 

tp = 1/8, p(z) = (z - 1)(1 + z 3 (l + z + z 2 )). 
Spectrum 

A » {0, 0.1550, 0.3913, 0.7037} + Z 

Coefficients 

o = 1 
ao.i550 ~ 1.55016 + 1.45286? 
ao.3913 « -0.505763 - 0.219617* 
ao.7037 « -0.1694 + 0.391761* 
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p(z) = (z-l)(l + z 3 (l + z)) A 

Figure 5.12. 
6. Questions 

Motivated by connections to problems from mathematical physics, we include a 
brief list of questions arising naturally in connection with properties of spectrum(PB) 
for the selfadjoint operators Pg considered here. Even though [Rob71] deals with 
second order operators, there are a number of parallels, for example both in [Rob71] 
and here, we are dealing with deficiency indices (2,2). And in both cases, we get 
that, when the boundary points are fixed, then the individual eigenvalues X n (B) 
from the list spectrum(Pe), for each n, depend continuously on B. 

(1) Robinson [Rob71]: one interval, but a second order operator, indices (2, 2), 
and a formula where the spectrum is a sequence of points in A (on the real 
line) arising by graph intersections. Our case: two intervals, and a first 
order operator, and the same conclusion as in Robinson, only with graphs 
of a different pair of functions. Open question: What can be said about 
"similar" of extensions of Hermitian operators with indices (2,2)? 

(2) Which self-adjoint restrictions of d 2 /dx 2 are not P 2 for some self-adjoint 
restriction P of id/dx? 

(3) Obtain a formula for the dynamics, i.e., description of t — > e ltPn f(x)? 

(4) Corollary 3.17 and plots of A„ as a function of n, suggests it might be 
interesting to obtain a quantitative estimate, depending on the other pa- 
rameters, on (5 = 5 W ? 

(5) More than two intervals. Some results in the paper do carry over to sets 
formed as union of n intervals, but not our detailed classifications. With n 
intervals, the preliminaries in Sections 1-2 still hold: the deficiency indices 
will be (n,n). The boundary form in Lemma 1.1 will also have a matrix 
representation. Corollary 2.3 will still hold; but instead there will be n left 
interval endpoints, and a corresponding set of right-hand side interval end- 
points. As a result, the boundary matrix B will now be an n by n complex 
unitary matrix; every B in the group U (n) will define a selfadjoint oper- 
ator Pg . The group U (n) allows an Iwasawa decomposition, and we will 
get an analogue of w in (2.10), only now it will be an positive selfadjoint 
matrix W, satisfying < W < I, referring to the usual order of selfadjoint 
matrices. One will also have a matrix version of the Mobius transform in 
(3.20) in Lemma 3.12. There is a theory of matrix operations, making use 
of fractional linear matrix operations. There is even a matrix version of the 
Poisson kernel (3.24). But our detailed classification in Theorem 3.14 will 
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not carry over. Nonetheless, each selfadjoint operator Pb will still have a 
pure point spectrum, but now with a much more subtle multiplicity config- 
uration. And, for the general multiple-interval case, it would be difficult to 
get a classification of all the spectral pairs anything close to what we have 
in our theorems in Section 5 above. 
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